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Abstract 

We introduce a generalized similarity analysis which grants a qualitative descrip- 
tion of the localised solutions of any nonlinear differential equation. This procedure 
provides relations between amplitude, width, and velocity of the solutions, and it is 
shown to be useful in analysing nonlinear structures like solitons, dublets, triplets, 
compact supported solitons and other patterns. We also introduce kink-antikink 
compact solutions for a nonlinear-nonlinear dispersion equation, and we construct 
a basis of finite wavelength functions having self-similar properties. 



PACS numbers: 05.45.-a, 47.20.Ky, 02.30.Px, ll.10.Lm 



1 Introduction 

Nonlinear dynamics, which is a relatively new field of study, is a very important frontier 
for probing natural phenomena. Active research efforts that focus on the mathematics 
and physics of nonlinear dynamical systems have emerged worldwide in various fields, 
including fluid dynamics, plasma physics, astrophysics, and even string theory Q. Among 
its unique features is the ability to describe a variety of patterns and particle-like 
traveling solutions ||. Other notable features of the theory include its description of 
solitons and breather modes, features in quantum optics [Q, molecular and solid state 
physics phenomena ||, and solitons in nuclear and particle physics H. 
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In contrast with linear theories which exhibit smooth regular motion, nonlinear mod- 
els require nonlinear partial differential equations (NPDE) and show strong couplings 
between different mechanisms and parts of the system. Also, the nonlinear interactions 
involve multiple scales [0] and are related with self-similar patterns or fractals. The NPDE 
solutions of physical interest are mostly localized and demonstrate good stability in time 
and through scattering with each other. Their shapes are related to the velocity, thus 
making the nonlinear patterns distinct from linear waves. In the asymptotic domain 
these solutions consist of isolated traveling pulses that are free of interactions. Close to 
the scattering domain, the nonlinear solutions obey nonlinear superposition principles. 

The main challenge in any nonlinear analysis is the construction of localized or finitely 
supported analytical solutions for the NPDE of interest. This challenge includes issues 
regarding the inexistence of a superposition principle for such solutions. Recent examples 
show that the traditional nonlinear tools (inverse scattering, group symmetry, functional 
transforms) are not always applicable Pf. On the other hand, from an experimental 
point of view one knows that patterns which are observed in nature - either stationary, 
growing, or propagating - generally have finite space-time extension and a multi-scale 
structure. Since soliton, even when they are localized, have an infinite extent, one needs 
other appropriate structures, and eventually self-similar bases. 

In this paper, wavelet-inspired approaches for localized solutions of NPDE are ex- 
plored. We propose a new similarity formalism for the qualitative analysis, and clasiffica- 
tion, of soliton solutions of nonlinear equations. This method provides relations between 
the characteristics of such solutions (amplitude, width and velocity) without the need 
of solving the corresponding NPDE. The method uses the multi-resolution analysis ||10| 
where traditional tools like the Fourier integrals or linear harmonic analysis are inadequate 
for describing the system. Wavelets are functions that have a space-dependent scale which 
renders them an invaluable tool for analyzing multi-scale phenomena. Wavelets have been 
used in signal processing, in problems involving singular potentials, pattern recognition, 



image compression, turbulence and even in radar and acoustic problems |L0[ [TTJ • More- 
over, the introduction of wavelet analysis in a study of NPDE is very natural because 
they can acommodate everything from strong variations, even singularities, to a smooth 
behavior. 

The second purpose of the paper is to show an example of construction of a nonlinear 
basis for a NPDE with nonlinear dispersion. There are many physical reasons favoring 
wavelets in the construction of nonlinear bases. For example, the breakup process of fluid 
drops has been shown to be self-similar, and in particular, the singularities (necks) look 
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identical at any scale [12[]. Other nonlinear oscillations of liquid drops, shells, bubbles, 
or even neutron stars involve such types of behavior 0, |9|. We introduce in section 2 
a qualitative similarity analysis for NPDE that yields relations between the amplitude, 
width and velocity of their traveling solutions, with many examples and predictions. In 
section 3, we introduce a nonlinear wavelet-like frame, associated with localized analytic 
solutions of a modified KdV equation, with nonlinear dispresion. 



2 Qualitative similarity analysis 

Finding analytic solutions for the nonlinear equations which describe physical phenomena 
is rather an exception than the rule. The traditional nonlinear tools like the inverse 
scattering theory or functional transforms are not always applicable || [16|]. Also, when 



phenomena of interest have many space (or time) scales, hundred times less or smaller 
than the dimension (or time scale) of the whole system, numerical methods may fail. 
This is the case of sharp propagating perturbations developing discontinuities, or shock 
waves. A simple option is the expansion of solutions in a basis of appropriate chosen 
linear modes. The most common choice are the Fourier series which have the advantage 
of orthogonality, but can hardly discriminate local behavior of phenomena. Moreover, 
information about the order of magnitude of the Fourier coefficients of a signal or wave 
u(x, t) is not sufficient for making conclusions about the size or scale of u. 

In these special situations, the use of bases of multiresolution analysis and wavelets 
has become popular ||13|| . We introduce a qualitative analysis for the localized traveling 



solutions belonging to any type of NPDE, in terms of Morlet continuous wavelet approach. 
In our case of interest, traveling localized perturbations, the most important information 
are provided by the scale of the pulse (or the width) denoted L, the amplitude A, and 
the group velocity V. The qualitative analysis in this section provides simple relations 
between these three parameters, without actually solving the equation. It also gives an 
estimation about the specific scale of the solutions. 

The procedure consistsin the substitution of all the following terms in the NPDE, 
according to the rule 

u t -> ±Vu x , u -> ±A, u x ->±A/L, u xx -> ±A/L 2 . . . , (1) 

and so forth for higher order of derivatives. Consequently, the NPDE is mapped into an 
algebraic equation in A, L and V. In Table 1 we present several examples of application 
of this substitution to some well-known and widely used NPDE in physics. 
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The validity of the method follows from the expansion of the soliton-like solution u(x) 



in Morlet wavelets [1C 



.2 



■$ a {x) =7r~ 1 / 4 e- iax -^, (2) 

where a describes the scale of this mother wavelet. The support of any Morlet wavelet 
is mainly confined in the (—1, 1) interval. We have the discrete Morlet wavelet expansion 
of u 

j k 

in terms of integer translations and dyadic dilations of the mother wavelet. In order 
to reduce the number of scales needed, the range of the summation should be choosen 
with an eye to the underlying physics. We use in the following the asymptotic formula 
describing the pointwise behavior of the Morlet wavelet series around of a point xq of 



interest fL4| . For a chosen xq and scale j, there is only one k and |e| < 1 such that the 
support of the corresponding ^ a ,j,k contains this point, k = 2 J x + e. We can express the 
solution and its derivatives in a neighborhood of this point 

u(x ) « *(-e) C j&xo+t = XX'( x o), 

3 3 

— {xq) « -i*(-e) £ ^aC j:2JxQ+€ = -ia ^ 2?Uj(x ), (3) 

for a chosen enough large compared to e. Since the coefficient l/a2 J represents the scale 
for each ^ a j,k Morlet wavelet, we can define it as a characteristic half-width Lj. And we 
finaly have in xq, from eqs.(3) 

where Uj(xo) = ^ a (— t)Cj,2ix +e ~ ^a(Q)Cj,vxo- Eq.(4) is the many scales generalization 
of the simpler formula in eq(l). With eq.(4) in hand we can investigate the structure of 
hypothetic soliton solutions, by choosing xq in the neighborhood of the maximum value of 
the solution, u(xq, 0) = A. Around this maximum, such solutions can be described very 
well by a unique scale L, and hence the solution and its derivatives can be approximated 
with the corresponding dominant term, by the substitutions in eqs.(l). 
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In Table 1 we present a series of examples of NPDE, identified in the first column by 
the name and the form of the equation. In the second column, we write the corresponding 
traveling localized solution, if an analytical form is known. Such solutions provide special 
relations between L, A and V, which are given in the third column. In the last column we 
introduce for comparison the results of this Morlet qualitative analysis, that is the relations 
between the three parameters, provided by eqs.(l). The usefulness of the approach may 
be checked, by a quick comparison between the fourth and the fith columns. While the 
results in the fourth column are possible only when one knows the analytical solutions, the 
results presented in the last column, obtained by the similarity approach, result directly 
from the NPDE, without actually solving it. 

The first line in the Table presents a linear case, for comparative purposes. The above 
Morlet wavelet approximation provides a correct expression for the dispersion relation 
(V = c — > k 2 = u 2 /c 2 ) with no constraint on either the amplitude A or on the width L. 

The case of the KdV equation is described in the second row of the Table 1. The 
method gives a general expression for L = L(A, V). For L to be related to A only, from 
the fourth column it results that the velocity V must be proportional to A. In this case 
we obtain exactly the well-known relation (column three) among the parameters in the 
exact solution. A prediction of the method is if we allow V to depend on a power of A. 
This means solutions with a higher nonlinear coupling between the shape and kinematics. 
A side effect would be a lower limit for A. Smaller solitons than this limit can move vith 
velocity proportional to the amplitude only. 

The same result is obtained for the MKdV equation (third row), except that in this 
case V needs to be proportional to the square of A in order to have L a simple function of 
A only. This prediction is again identical with that in the exact solution (third column). 
Moreover, the same relations remain valid even for the new exotic solutions of the MKdV 



equation of compacton type, [15] 



u(x, t) 



32k cos[k(x - AkH)} 2 
3(1- § cos[fc(z - 4fc 2 t)] 2 ) ' 



which has L = 5n/6k or 7r/6fc, that is L ~ 1/A like in the Table 1. 

Next example (4-th row) is provided by a generalised KdV equation, in which the 
dispersion term is quadratic 

Vt + {i] 2 ) x + ( V 2 ) xxx = 0. (5) 



5 



Eq.(5), known as K(2,2) from the two quadratic terms, admits compact supported trav- 
eling solutions, named compactons [§, [16], [17], |I8| . The compactons are powers of trigono- 
metric functions defined on a half-period, and zero otherwise. In general, they have the 
form Acos a d(x — ct), and different from solitons, their width is independent of the ampli- 
tude. This is the fact that provides a connection with wavelet bases. They are characterize 
by a unique scale, and it is this feature that makes it possible to introduce a nonlinear 
basis starting from this "mother" function. For eq.(5) the compacton solution is given by 



4V 

7] C (X - Vt) = ~Y C0 ^ 



X 



Vt 



(6) 



if \x — Vt\ < 2tt and zero otherwise. Here the velocity is a function of the amplitude. 
Notice that the width L = 4 of the wave is independent of the amplitude. The quadratic 
dispersion term is characteristic for the nonlinear coupling in a chain. 

The general compacton solution for eq.(5) is actually a "dilated" version of eq.(6). 
That is, a combination of the first rising half of the squared cos in eq.(6), followed by a 
flat domain of arbitrary length (A), and finaly followed by the second, descending part 
of eq.(6). Actually, this combination is just a kink compacton joined smoothly with an 
antikink one 



Vkak(x - Vt] A) 



COS 



0... 

4V 
3 

4V 

3 



f cos 2 



-vt 



-2tt < x - Vt < 



< x - Vt < A 



(7) 



-Vt-X 



X<x-Vt<X + 2iT 



0... 



In Fig. 1 we present a compacton, eq.(6), a kink-antikink pair (KAK) described by 
eq.(7), both with the same amplitude and velocity. Although the second derivative of this 
generalized compacton is discontinuous at its edges, the KAK, eq.(7), is still a solution 
of eq.(5) because the third derivative acts on u 2 , which is a function of class C3. Finally, 
we can construct solutions by placing a compacton on the top of a KAK, like in the third 
solution in Fig. 1. Such a solution exists only for a short interval of time, since the two 
structures have different velocities. The solution is given by 



rj(x,t) = 7] kak (x - Vt; A) + [7] c (x - Vt 



4V\ 
2vr) + ^)x( 



x 



V't-2n, 



2vr 
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for < t < (A — 47r) / (V — V) . Here x( x ) is the support function, equal with 1 for \x\ < 1 
and in the rest, and V = 3max{r] c }/4 + 2V. 

For the K(2,2) compacton, eq.(6), the exact relations between the parameters are 



A = AV/3 and L = 4, fl7| . The relation provided by the similarity method, in the 
last column of the forth row, predicts the existence of the compacton. That is, for a 
linear dependence between the amplitude and the speed, the half-width is constant and 
does not depend on A or V. This fact (L = L = const.) is a typical feature of K(2,2) 
compactons. Moreover, in literature there was found numerically that for any compact 
supported initial data, widder than Lq, the solution decomposes in time into a series of Lq 
compactons, Fig. 2. For narrower initial data the numeric solution blows up. There is no 
exact or analytic explanation of this effect, so far. The similarity method can give a hint 
in this situation, too, by using the graphic of the relation L = L(V, A) provided by this 
qualitative method. In Fig. 3 L is ploted versus V, for several values of A (larger values 
of A translate the curves to the right). The half- width of a stable compacton was chosen 
Lq = 0.707. Above this value, Fig. 3a, wider compact pulses produce an intersection for 
each curve (each A) with the axis Lq providing series of compactons of different heights, 
like in the numerical experiments, Q, |T7|]. Below this Lq line, all the curves approache 
infinite amplitude, providing instability of narrower shapes. 

Another good example of prediction of the method is exemplified in the case of a 
general convection- nonlinear dispersion equations, denoted K(n,m) 

Vt + (V% + = 0. (9) 

Compacton solution for any n ^ m are not known in general, except some particular 
cases. In this case we find a general relation among the parametrs, for any n, m, shown 
in the 5 th and 6 th rows. These general relations L(A, V) approache the known relations 
for the exact solutions, in the particular cases like n = m (5 th row), n = m = 2 (4 th row), 
n = m = 3. And n = 3, m = 2; n = 2, m = 3 in the 6 th row. These results can be used to 
predict the behavior of solutions for all values of n, m. 

Similar analysis can be done in the case of sine-Gordon equation, if we ask that velocity 
be proportional with L 2 (7 th row). In this case we obtaining a transcendental equation in 
A, which is just the case of the sine-Gordon soliton. In the 8 th row, we present the cubic 
nonlinear Schrodinger equation (NLS) which has a soliton solution, too M. This equation 
arises, for example, in nonlinear optics or in the polaron model in solid state physics, 
[[|, In the general case of a NLS of order n (9 th row), when the general analytical 
solution is unknown, the method predicts a special L = L(A, V) dependence shown in 
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the fourth column and in Fig. 4. Contrary to third order NLS, where the dependence 
of L with A is monotonous for V =~ ±A (n = 3 in Fig. 4), at higher order, the L(A) 
function has a discontinuity in the first derivative. This wigle of the function (Fig. 4 for 
n = 4) yields at a critical width, producing bifurcations in the solutions and scales. As 
a consequence, initial data close to this width can split into doublet (or even triplet for 
higher order NLS) solutions, with different amplitudes. Such phenomena have been put 



into evidence in several numerical experiments for quintic nonlinear equations |L7|, 

In the following, we present another example of applications of this qualitative ap- 
proach, related to a new type of behavior of nonlinear systems. Traditional solitons move 
with constant speed on a rectilinear path (except for the roton, M which has a circular 
trajectory with constant angular velocity). The speed is usually equal with the amplitude 
scaled with a constant. Higher solitons travel faster and there are no solitons at rest (zero 
speed asks for zero amplitude). They can travel in both directions with oposite signs 
for the amplitude. The situation is different in the case of compactons, which allow also 
stationary solutions. When linear and nonlinear disspersion occur simultaneously, like in 
the so called K(2,l,2) equation 

u t + (u 2 ) x + (u) xxx + e{u 2 ) xxx = 0, 

where e is a control parameter, the similarity approach yields a dependence of the form 



L = y/(±A + €)/(V±A), 

which still provides a constant width if V — ±A+2e. In this case the speed is proportional 
with the amplitude, but can change its sign even at non-zero amplitude. Solutions with 
larger amplitude than a critical one (A cr i t = =F2e) move to the right, solutions having the 
critical amplitude are at rest, and solutions smaller than the critical amplitude move to 
the left. This behavior was explored in fPTf . However, such a switching of the speed is 
not necessarily a feature of the nonlinear dispersion. A compacton of amplitude A on the 
top of a KAK solution of amplitude S 

u(x, t) = Acos 2 ( - ~ A Vt ) + 5, (10) 



is still a solution of the K(2,2) equation, u t + {u 2 ) x + {u 2 ) xxx = 0, with the velocity given 
by V = 1 1 28 + A ) . For A = —25 the compacton becomes a stationary anticompacton, 
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embedded in the moving, supporting KAK. Such an example is presented in Fig. 5 
for a slow-scale time- dependent amplitude compacton. The induced oscillations in the 
amplitude transform into oscillations in the velocity. While not the topic of this paper, 
such a dynamic system has been analysed and it will be published soon. The key to 
such a conversion of oscillations is the coupling between the traditional nonlinear picture 
(convection-dispersion-diffusion) and the typical Schrodinger terms. 

A last application of this method, occurs if the KdV equation has an additional term 
depending on the square of the curvature 

u t + uu x + u xxx + e(u 2 xx ) x = 0. (11) 

This is the case for extremely sharp surfaces (surface waves in solids or granular materials) 
when the hydrodynamic surface pressure cannot be linearized in curvature. Such a new 
term yields a new type of localized solution fulfilling the relations 



4eA 



y\ ±y/l -SeA(A±V) - 1 

If we look for a constant half- width solution (compacton of 1/L = a) we need a dependence 
of velocity of the form V = (l + a 2 e/8)A + l/8eA + a/4:. There are many new effects in this 
situation. The non-monoton dependence of the speed on A introduces again bifurcations 
of a unique pulse in dublets and triplets. Also, there is a upper bound for the amplitude 
at some critical values of the width. Pulses narrower than this critical width drop to zero. 
Such bumps can exist in pairs of identical amplitude at different widths. They may be 
related with the recent observed "oscillations" in granular materials, |T7|, [I9fl . 

As the examples presented in Table 1 proved the above method provides a reliable cri- 
terium for finding compact suported solutions. The reason this simple prescription works 
in so many cases follows from the advantages of wavelet analysis on localized solutions. 
We stress that this method has little to do with the traditional similarity (dimensional) 
analysis |8], [12|, |1(| [17|, [18, |2(J. In the latter case one obtains relations among powers of 



A, L and V, not relations with numeric coefficients like those found in our method. 
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3 Compacton kink-antikink pairs and the multireso- 
lution frame 



A common feature of all NPDE and of the finite differences equations is the existence 
of compact supported solutions. Compactons and discrete wavelets are typical examples. 
An interesting general conclusion can be obtained if we look at a one-dimensional model 
described by the most general NPDE dynamical equation 

d t u = 0{x, d x )u, (12) 

where O is a nonlinear differential operator. By taking into account only traveling so- 
lutions, this NPDE reduces to a NODE in the coordinate £ = x — Vt for an arbitrary 
velocity V. If is a compact solution it results that it is not unique for given initial 
compact data. If one chooses zero initial value for the solution and its derivatives up 
to the requested, in a certain point £ of the £ axis, these conditions can be fulfilled by 
any linear combination of disjoint translated versions of one particular solution, placed 
everywhere on the axis except £o- Consequently, for such initial data, the solution is not 
unique. This result shows that the compact supported property of the initial data and 
the solution, implies its non-uniqueness. 

Since we can transform the NODE into a nonlinear differential system of order one 

^ = F(£,U), U=(u,d x u,...), (13) 

we can apply the fundamental theorem of existence and uniqueness to solutions of eq.(13), 
for given initial data ?7(£o) = Uq. If the function F in eq.(13) fulfills the Lipschitz 
condition (its relative variation is bounded) than, for any initial condition, the solution is 
unique, pl| . Since any linear function is analytic and hence Lipschitz, we conclude that 
only nonlinear functions F allow the existence of compact supported solutions. Thus, 
a compact soliton implies non-uniqueness in the underlying NPDE, which implies non- 
Lipschitzian structure of the NPDE and hence the existence of nonlinear terms. 

In the following we investigate some compact solutions of the K(2,2) equation. The 
high stability against scattering of the K(2,2) compactons, or compacton generation from 
compact initial data, suggest they may play the role of a nonlinear local basis. We know 
form many numerical experiments, || |TB, |T7|, |T8|] , that any positive compact initial data 



decomposes into finite series of compactons and anticompactons. This suggests that an 
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intrinsic ingredient for a nonlinear basis could be the multiresolution structure of the 
solutions, similar with the structure of scaling functions in wavelet theory. 

The compactons given in eqs.(6,7) have constant half- width and hence describes a 
unique scale, which can cover all the space by integer translations. From the point of 
view of multi-resolution analysis, the K(2,2) equations acts like a L-band filter, allowing 
only a particular scale to emerge for any given set of initial condition. To each scale, from 
zero to infinity, we can associate a K(2,2) equation with different coefficients. However, 
the compacton solution is not the unique one with this property. For a given K(2,2) 
equation, we can thus extend the scale from L to any larger scale. These more general 
compact supported solutions are still C^R) and are combinations of piece- wise constant 
and piece-wise cos 2 functions. The simplest shape is given by a half-compacton prolonged 
with a constant level, that is a kink solution. The basis solution is a kink-antikink (KAK) 
compact supported combination, Fig. 1. Such kink-antikink pairs of different length, can 
be associated with other compactons, or KAK pairs, one on the top of the other 



Vcomp+KAK(x ~ Vt] A) 



0... 



v-Vt 



f cos 2 
f, 0<, 



-2tt < x - Vt < 



Vt<5 



-Vt 



f + l(V'-2V) cos 

5 + Air < x — Vt < A 



, 5<x-Vt<5 + 4n (14) 



f cos 2 
0... 



v-Vt-X 



\<x-Vt<\ + 2n 



where 5 < A characterizes the initial position (at t — 0) of the top compacton, with respect 
to the flat part of the KAK solution. The amplitude 4(V — 2V)/3 of the compacton, and 
the amplitude 4V/3 of the KAK, are related to their velocities V and V, respectively. 
The length of the flat part, A, is arbitrary. A compound solution is not stable in time since 
the different elements travel with different velocities. The total height of the compacton 
is 4(V — V)/3. Since the higher the amplitude is, the faster the structure travels, the 
top compacton moves faster than the KAK, and at a certain moment it passes the KAK. 
Because the area of the solution is conserving, such a compound structure decomposes 
into compactons and KAK pairs. Similar and even more complicated constructions can 
be imagined, with indefinite number of compactons and KAK's, if one just fulfills the C3 
continuity condition for the square of the total structure. Such structures, defined at the 
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initial moment can interpolate any function, playing a similar role with wavelets or spline 
bases. It has been also proved that the KAK solutions are stable, by using both a linear 
stability analysis and Lyapunov stability criteria, (IB], |2D 



For a given K(2,2) equation, the compacton solution, eq.(6) and in addition the family 
of KAK solutions, eq.(7) can be organized as a scaling functions system. They act like a 
low-pass filter in terms of space-time scales and give the opportunity to construct frames 



of functions from the wavelet model, [10, 11, [L 



For the sake of simplicity we will renormalize the coefficients of the K(2,2) equa- 
tion such that the support of the simple compacton is one. That is, we take r) c (x,t) = 
Vkaki^ix — Vt), 0) on the interval \x — Vt\ in [—1/2, 1/2]. We construct a multiresolution 
approximation of L 2 (R), that is an increasing sequence of closed subspaces Vj, j G Z, of 
L 2 (R) with the following properties, |TT| , |I3"| 

1. The Vj subspaces are all disjoint and their union is dense in L 2 (R). 

2. For any function / G L 2 (R) and for any integer j we have f(x) G Vj if and only if 
D^ 1 f{x) G Vj_i where D~ l is an operator that will be defined later. 

3. For any function / G L 2 (R) and for any integer k, we have f(x) G V is equivalent 
with f(x - k) G V . 

4. There is a function g(x) G Vq such that the sequence g(x — k) with k G Z is a Riesz 
basis of Vq. 

In the case of compact solutions of K(2,2) of unit length, we chose for the space Vq that 
which is generated by all translation of r] c with any integer k. The subspaces Vj for 
j > are generated by all integer translations of the compressed version of this function, 
namely, by r/kaki^^ix — Vt), 0). The subspaces Vj for j < are generated by all integer 
translations of the KAK solution of length A2 J — 1. For example, VI i is generated by 
Vkaki^ix — 2 J Vt),0). The spaces Vj,j > are all solutions of K(2,2); the others are not. 
The function g(x) is taken to be rjkak{^(x — Vt), 0). It is not difficult to prove that these 
definitions fulfill restrictions one, three, and four. As for the second criterion, we define 
the action of the operator D^f{x) = f(2x) if f(x) G Vj with a j positive integer, and 
D-'rikakiTrVix - VVt),2-i - 1) = %afc (7r2^'(x - 2^ +1 Vt),2^ +1 - 1) for negative j. In 
conclusion, we construct a frame of functions made of contractions of compactons and 
sequences of KAK solutions. We can write the corresponding two-scale equation which 
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connects the subspaces (the equivalent of eq.(15)), 

Vkak(n(x - Vt), 1) = r] kak (n(x - Vt), 0) + r] kak (rr(x - Vt - 1), 0). (15) 

We will denote generically by r] k j the elements of this frame, that is 

Vk,j( x ) = Vkak(n(x - 2 3 Vt - k), 2 3 - l)|t= , 

where t — means that we neglect the time evolution, but the amplitude is still amplified 
with a factor of 2 J , in virtue of relation r] max = 4V/3. In the following, we can expand 
any initial data for the K(2,2) equation in this basis. 

M x ) = J2J2 CkjVkAx). (i6) 

k j 

We notice that the following equality holds for j',j 

f -t0 k' = k- 2?'-i, .., (k + 1) • 2?'-i - 1 n7 , 
""M*' { = otherwise. (17) 

After some rather elaborate algebraic calculations and by using eq.(21), we show that the 
square of this function (since the equations is nonlinear and of order two) will be given 
by 

« 2 ( i ) = EEE CkjCk'j' 

k,j j'>j k'el 

,111 , 

X EE- E Va(i 1 ,i 2 ,...,i j ,_ j ), j ')Vk',j', (18) 

v n=oi 2 =o V-j =0 

where / is the range of k! described in the first line of eq. (21), and 

a(t 1 ,t 2 ,...,t j/ _ j ) = Y^iil 3 3l+ 2 



i 



+k2 u '~ j)j 

From eq.(21) we notice that in eq.(22)the unique nonzero terms are those for which 
a(ii, iii ■■■■> ij'-j) = k' with k! e /. This result express the following simple fact. The 
initial data is expanded in different scales and different translations. The translations are 
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mutually orthogonal so they do not give a contribution to the square. When we have to 
multiply two different scales in the expression of the square, we reduce the wider scale in 
terms of linear combination of the narrower one by using the two-scale equation, eq.(19). 
This is what eq.(22) expresses. Out of all the terms in such a product only approximately 

give non-zero contributions. In other words this number is 
given by the number of solutions of equation a(i±, 12, ij'-j) = k', with k! e /. This 
is the primary advantage of treating nonlinear problems with a basis that has a scale 
criterion. Another advantage is that all the function in the basis are actually contractions 
or dilations, and translations of only two basic ones. 

4 Comments and conclusions 

In the present paper we introduce new applications for wavelets, in the field of the study 
of localized solutions of nonlinear differential equations. The existence of compactons and 
discrete wavelets underlines a common feature of NPDE and finite differences equations, 
that is the existence of compact supported solutions. We propose a new similarity formal- 
ism for the qualitative analysis and clasiffication of soliton solutions, without the need of 
solving the corresponding NPDE. Also, we proved that starting from any unique soliton 
solution of a NPDE, we can construct a frame of solutions organized under a multires- 
olution criterium. This approach provides the possibility of constructing nonlinear basis 
for NPDE. We show that frames of self-similar functions are related with solitons with 
compact support. In addition, we notice the evidence that compactons fulfil both char- 
acteristics of solitons and wavelets, suggesting possible new applications. Such unifying 
direction between nonlinearity and self-similarity, can bring new applications of wavelets 
in cluster formation, at any scale, from supernovae through fluid dynamics to atomic and 
nuclear systems. The similarity approach can be applied with succes to the physics of 
droplets, bubbles, traveling patterns, fragmentation, fission and inertial fusion. 

Supported by the U.S. National Science Foundation through a regular grant, No. 
9970769, and a Cooperative Agreement, No. EPS-9720652, that includes matching from 
the Louisiana Board of Regents Support Fund. 
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Figure Captions 



Fig. 1 

A compacton and a kink-antikink pair solution (KAK) of the equation K(2,2), both 
having the same amplitude, and hence velocity V. To the right, there is a smaller 
compacton on the top of KAK. The upper compacton has higher speed, V. 

Fig. 2 

A finite series of K(2,2) compactons emerging from initial compact data, with the 
width larger than the compacton width. 

Figs. 3 

The half-width L versus velocity V for the K(2,2) equation, for different amplitudes 
A. Figure 3a shows widths larger than L compacton = 3/4, and Fig. 3b shows narrower 
widths, L < 3/4. Amplitude increases from left to right, in the range 0.01-0.85. 

Fig. 4 

The half-width L plot versus amplitude V, for the third (n=3) and forth (n=4) 
order NLS equation, in two V = ±A cases. From the figure one notes that the 
quartic NLS equation yields bifurcations in the solutions. 

Fig. 5 

The solution of the mixed linear plus nonlinear-dispersion K(2,2) equation, in the 
case of a solution with a slow oscillating shape. 
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Table 1: Nonlinear equations, exact solutions and Morlet similarity analysis. 



Equation 



Solution 



Relations 



Wavelet 



Linear wave 
u xx - {l/c 2 )u tt = 



£ C k e^ kx± ^ k 2 = u 2 /c 2 V = c 

A,L arbitrary 



KdV=K(2,l) 
u t + 6uu x + u xxx = 



A sech 2 ^P 



T - 2 
L ~ V A 

V = 2A 



y/\±V±6A\ 



MKdV=K(3,l) 

u t + u 2 u x + u xxx = 



A sech^P 



L = 1/A 
A = W 



^\±V±6A 2 \ 



K(2,2) 
u t + {u 2 ) x + (u 2 ) 



A cos 



2 x-Vt 



L = 4 



8A 



if 



L ^ -c ^ y |±V±2A| • 

|(x-V^)/L| < tt/2 V = 3,4/4 V = -3A/2, L = 4 



K(n,n) 

«t + + (« n )xxx = 



2 / 





A = 


2Vn 




n+1 


2nn 


L = 


An 


n-1 


(n-1) 



and else 



t(n 2 +l) 



±a±n 

if V = aA^ 1 



K(n,m) unknown 
u t + (u n ) x + (u m ) xxx = iin^m 



n(n 2 +l)A n - 
±V±mA m - 



sine-Gordon 
u xt — sin u = 



A tan 7 e £ 



A = 4 



±jj = sin A 



NLS(3) 
% + V xx + ^ 3 = 



z-vt 



^ = i 
A ~ V 



±V±^\V 2 -4A? 
2A 2 



NLS(n) 
% + ^ xx + V n = 



unknown 



L = 



±v± v /|y 2 -4A™| 

2A™ 



